Dynamical and sequential decay effects on isoscaling and density dependence of the 

symmetry energy 
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The isoscaling properties of the primary and final products are studied via isospin dependent 
quantum molecular dynamics (IQMD) model and the followed sequential decay model GEMINI, 
respectively. It is found that the isoscaling parameters a of both primary and final products keep 
no significant change for light fragments, but increases with the mass for intermediate and heavy 
products. The dynamical effects on isoscaling are exhibited by that a value decreases a little with 
the evolution time of the system, and opposite trend for the heavy products. The secondary decay 
effects on isoscaling are reflected in the increasing of the a value for the final products which 
experiences secondary decay process. Furthermore the density dependence of the symmetry energy 
has also been explored, it is observed that in the low densities the symmetry energy coefficient 
has the form of C sy m(p) ~ Co(p/po) 7 , where 7 = 0.7 ~ 1.3 for both primary and final products, 
but Co have different values for primary and final products. It is also suggested that it might 
be more reasonable to describe the density dependence of the symmetry energy coefficient by the 
C sym {p/po) ~ Ci(p/po) 7so/t + C 2 (p/po) 7stl// with y so ft < 1, Jstiff > 1 and Ci,C 2 constant 
parameters. 

PACS numbers: 24.10.Nz, 24.10.-i, 25.70.Pq, 25.70.-z, 24.80.+y 



I. INTRODUCTION 

The isotopic composition of the nuclear reaction prod- 
ucts Q contains important information on the role of 
the isospin on the reaction process. Recently increased 
interest in the N/Z degree of freedom and its equilibra- 
tion, as well as the isospin asymmetry dependent terms 
of the nuclear equation of state (EOS) |E IE EL IE IE ' 
has motivated detailed measurements of the isotopic dis- 
tributions of reaction products. It has been shown that 
isospin effects can be studied by comparing the yields of 
fragments from two similar reactions that differ only in 
the isospin asymmetry jE0> m * ms case i the effect of 
sequential decay of primary fragments can be bypassed 
to a large extent. It has been revealed that for statistical 
fragment production mechanism(s), if two reactions oc- 
curring at the similar temperature have different isospin 
asymmetry, the ratio R.2i(N, Z) of the yields of a given 
fragment N and Z obtained from the two reactions 2 and 
1 exhibits an exponential dependence on N and Z of the 
form [lEEEI 

R 21 (N,Z)=Y 2 (N,Z)/Y 1 (N, Z) = C exp{aN + 0Z), 

. (1) 

where a and (3 are two scaling parameters and C is an 
overall normalization constant. This behavior is called 
isoscaling [TJj. Isoscaling has been obtained in different 
reactions and theoretical calculations, such as in evapora- 
tion reactions [l^, deep inelastic reactions [T3IIT5I ]. fission 
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reactions 0] and multifragmentations [E 0, 0, 113 , 
where is the biggest amount of experimental isoscaling 
data comes from. Isoscaling phenomena in different re- 
action mechanisms have also been investigated by various 
theoretical models, such as Langevin equation combining 
with the statistical decay model for the fission dynam- 
ics 0, |2(]], isospin-dependent Lattice Gas model [2l| 
for the multifragmentations, Anti-symmetrized Molec- 
ular Dynamics (AMD) model [2E], I sosp in dependent 
Quantum Molecular Dynamics model [23, Microcanon- 
ical Statistical Multifragmentation model, Expanding 
Emitting Source model, Canonical model Micro- 
canonical Multifragmentation model J2J] and so on. 



The isospin dependence of the nuclear equation of state 
(EOS) is one of the most important properties in nu- 
clear matter and reactions, in particular in nuclear as- 
trophysics, supernova and neutron stars l25j |. Al- 
though the nuclear symmetry energy at normal nuclear 
matter density po = 0.16/to -3 has been determined 
to be around 30 MeV from the empirical liquid-drop 
mass formula [2(|, however, its values at sub- and super- 
densities are poorly known. Studies based on various 
theoretical models also give widely different predictions 
0, 0, 11E ■ Multifragmentation is generally consid- 
ered as a low density phenomenon with a high degree of 
thermalization which is believed to be reached. It has 
been shown [111 llq that the isoscaling parameter a is 
directly related to the coefficient C sym of the symmetry 
energy term of the nuclear binding energy, the follow- 
ing relation has been obtained both in the framework of 
the grand-canonical limit of the statistical multifragmen- 
tation model |Tl| and in the expanding-emitting source 
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where Z±, A\, and Z2, Ai refer to the charge number and 
mass number of the fragments from reactions 1 and 2 re- 
spectively. Using this relation, from the extracted values 
of a and T of the fragments in the reaction, symmetry 
energy coefficient C sym in the EOS could be derived. 

The paper aims to investigate isoscaling in multifrag- 
mentation reactions by the dynamical IQMD model and 
the afterburner sequential decay model GEMINI calcula- 
tion, explore the dynamical and secondary decay effects 
on the isoscaling phenomenon. Different isoscaling prop- 
erties between light and heavy products have revealed. 
Through the relationship of equation between isoscal- 
ing parameter a and symmetry energy coefficient C sym , 
the form of symmetry energy density dependence can be 
deduced. Section El gives a brief review of the IQMD 
and GEMINI model. Section EH presents the IQMD 
and IQMD+GEMINI calculation results on isoscaling ob- 
tained from the different system evolution time, primary 
and final products, discusses the properties of the isoscal- 
ing parameter a. Section llVl investigates the density de- 
pendence of the symmetry energy coefficient C sym for 
the primary and final products following some results in 
section Hill The conclusions are drawn in section Ivl 



states. Through a Wigner transformation of the wave 
function, defining L = of and cr 2 — h 2 /4a 2 = h 2 /4L the 
iV-body phase-space distribution function is given by a 
symmetry form in coordinate and momentum space 
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here i"i and pi are the center of the ith wave packet 
in the coordinate and momentum space, ay and a p are 
the widths of wave packets in coordinate and momen- 
tum space, respectively. The Wigner representation of 
the Gaussian wave packet obeys the uncertainty relation 
a r a p = ft/2. The densities in coordinate and in momen- 
tum space are given by 
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In the IQMD model, the nucleons in a system move 
under a self-consistently generated mean field, and the 
time evolution of and is governed by Hamiltonian 
equations of motion 



II. MODEL OVERVIEW 



A. Dynamical model: IQMD 
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The Quantum Molecular Dynamics (QMD) model ap- 
proach is an 71-body theory to describe heavy ion reac- 
tions from intermediate energy to several GeV/nucleon, 
it is classical in essence because the time evolution of the 
system is determined by classical canonical equation of 
motion, however many important quantum features are 
included in this prescription. At intermediate energies 
the heavy ion collision is mainly governed by three com- 
ponents: the mean field, two-body collision and Pauli 
blocking. General review about the QMD can be found 
in reference of Aichelin [2|j tog ether with several other 
similar models such as AMD [311 and Constrained Molec- 
ular Dynamics (CoMD) model |32|. 

The IQMD model is based on the general QMD to 
include explicitly isospin-degrees of freedom ■ In the 
QMD model, each nucleon is represented by a Gaussian 
wave packet, 
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where v\ is the center of the ith wave packet in the coor- 
dinate space, and L is the so-called Gaussian wave packet 
width (here L = 2.16/to 2 ). The total TV-body wave func- 
tion is assumed to be the direct product of these coherent 



The Hamiltonian H which consists of both kinetic en- 
ergy and effective interaction potential energy is given 
by 



H 
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where fi is the mass of a nucleon, U dd the density- 
dependent (Skyrme) potential, JJ yuk the Yukawa poten- 
tial, JJ sym the symmetry energy term and JJ coul the 
coulomb energy, they have the following forms 
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TABLE I: The parameters adopted in present work. 
a{GeV) a PjGeV) 7 pojfm' 3 ) v y (GeV) m C sym {GeV) 



-0.356 0.303 1.17 0.16 



-0.0024 0.83 



0.032 



"Here a, /3 and 7 values correspond to soft potential 
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Here = |r*i — r*j- 1 is the relative distance of nucleon 
i and j, Ti z the 2th component of the isospin degree of 
freedom for the ith nucleon, which is equal to 1 and -1 for 
proton and neutron, respectively, C sym is the symmetry 
potential coefficient. The parameters in this work are 
listed in Table [Q 

The nucleon-nucleon (NN) cross section used in 
present IQMD model is experimental parametrization 
|33j (<Jexp) which is isospin dependent, the neutron and 
proton are distinguished from each other in the initializa- 
tion of projectile and target nuclei, and the Pauli blocking 
of proton and neutron in IQMD is also treated separately 
as well 30] . In our work, the fragment (or cluster) is rec- 
ognized by a simple coalescence model: i.e., nucleons are 
considered to be part of a cluster in every moment at 
least another nucleon is closer than r m i n < 3.5 fm in the 
coordinate space. 



B. Sequential decay model: GEMINI 

The primary fragments from IQMD are excited, in or- 
der to compare with the experimental observable, the 
GEMINI sequential decay model |3J, |35| has been used 
as an afterburner to follow the de-excitation of these ex- 
cited fragments. In GEMINI all possible sequential decay 
channels from light particle evaporations to symmetry 
fission are considered. The sequential decay is traced by 
a Monte-Carlo technique until the excitation energy of 
the excited fragment is exhausted. In IQMD calculation 
the angular momentum is not considered, hence the sec- 
ondary decay treatment in the GEMINI does not include 
the angular momentum calculation. In other word, for 
each sequential decay tracing, the initial angular momen- 
tum is set to zero in GEMINI, and in present work all 
asymmetric divisions are considered. More details about 
GEMINI can be found in Refs. [HHU. IQMD calcula- 
tion is on the basis of event-by-event, so the secondary 
decay by GEMINI is managed event-by-event too. Each 
fragment calculated in one IQMD event will put into the 



GEMINI to let them decay. The excitation energy of the 
primary fragments is calculated by following equation 
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where Ef^}- is the incident kinetic energy in Center of 
Mass (cm.) system and E?^ is the fragment kinetic 
energy in cm. system. The summer is made over all the 
fragments of each event, the Q- value is calculated by the 
mass excess in each event, then the excitation energy per 
nucleon E* /A was uniformly distributed in all nucleons. 

In this paper all asymmetric divisions was consid- 
ered(including both intermediate mass fragment(IMF) 
and fission decay) in the GEMINI secondary decay calcu- 
lation, after GEMINI calculation, the final products were 
collected inclusively again. 



III. ISOCALING PROPERTIES AND THE 
DYNAMICAL AND SECONDARY DECAY 
EFFECTS 

Calculations were carried out in two similar central 
collisions with impact parameter b — 1 of 40 Ca + 40 Ca 
and 48 Ca + 48 Ca at incident energy E/A = 35 MeV. 
In order to investigate the isoscaling dynamical effect, at 
different reaction time fragments were recorded until 400 
fm/c, about 50 thousands events are simulated. 

Equation Q can be expressed as a function of only 
atomic number Z when the neutron number N is fixed 

R*i{Z) = Y 2 (N,Z)/Y 1 (N,Z) = C'exp(aiV), (14) 

or a function of neutron number N when the atomic num- 
ber Z is fixed 

Rn{Z) = Y 2 (N,Z)/Y 1 {N,Z) = c" exp(J3Z). (15) 

Isoscaling of hot fragments was observed at different 
reaction time from t = 160fm/c to 400fm/c, fragments 
mentioned here include all of the fragments formed in 
that moment. In fig. 1 two examples of R21 were shown 
for the fragments at reaction time t = 240 fm/c as a 
function of atomic number Z and neutron number N, 
respectively. R21 shows perfect linear dependence as a 
function of N or Z from light fragments to heavy frag- 
ments. Here only reaction time t = 240 fm/c are plotted 
as one sample, the results at other reaction times have 
similar behaviors. Here we adopt the widely used conven- 
tion to denote with the index " 2" the more neutron-rich 
system 48 Ca + 48 Ca and with the index "1" the more 
neutron-poor system 40 Ca + 40 Ca. 

From Equation (|14fl and (|15l) the isoscaling slope pa- 
rameters a and j3 can be directly extracted, in this paper 
only parameter a is analyzed, since parameter (3 has sim- 
ilar trend with a except that f3 value is negative. In fig. 
Pfa). parameter a as a function of the fragment charge 
number Z was printed at several reaction time, one can 
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FIG. 1: Isotopic yield ratios (upper panel) and Isotonic yields 
ratios (bottom panel) of primary fragments at reaction time 
t = 240 / m/c between 48 Ca + 48 Ca and 40 C 'a + 40 Ca in the 
IQMD model with impact parameter b = 1/m and incident 
energy Ei nc — 35MeV/A, different symbols from left to right 
represent the isotope sequence from Z — 1 to Z — 23 (upper 
panel) and isotone sequence from N — 1 to Z = 25 (bottom 
panel), in this figure the error bar comes from only statistics 



find that for the light and intermediate mass fragments, 
such as Z < 15, a value does not change too much with 
the fragment charge number Z, but for the heavy frag- 
ments, a shows different behavior from the light ones, it 
increases with the fragments charge number Z. One can 
also find that the isoscaling parameter a presents very 
small decrease when the reaction time becomes longer, 
at t = 140 fm/c a lies on the top, t = 400 fm/c a locates 
at the bottom. However, this is suitable only for the 
light products, and it is opposite for the heavy products 
though the dynamical effect is not clear. 

From a practical point of view, the isotopic distribu- 
tions of a given Z fragment Y (N) \z can be approximately 
described by Gaussian function. Using the general nota- 
tion for the isotopic distribution in a given Z, Y(N)\z 



FIG. 2: Isoscaling parameter a as a function of the fragment 
atomic number Z, different symbols represent a at different 
reaction time, panel (a): IQMD calculation results; panel (b): 
IQMD+GEMINI calculation results. 



can be described by a single Gaussian function as 



Y(N)\ Z = Cexp 



(N-(N z )f 
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where (Nz) is the centroid of isotopic distributions, and 
az describes the variance of distributions for each ele- 
ment of charge Z , then for a fixed element of charge Z 
the ratio of isotopic yields writes 
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where (Ni) and (A^) are the centroids for isotopic dis- 
tributions of neutron-poor system 40 Ca + 40 Ca and 
neutron-rich system 48 Ca + 48 Ca, a\ and a 2 describe 
the variance of distributions in these two systems, re- 
spectively, here the fragment charge number Z is fixed. 
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FIG. 3: The centroid (N) and width a 2 of Gaussian isotopic 
distribution as a function of the atomic number Z at different 
reaction times for the primary products in IQMD calculation 
and final products in IQMD+GEMINI calculation. Panels 
(a) -(f): centroid (N) as a function of Z at different evolution 
time in two systems; panel (e) and (f): a 2 as a function of Z, 
(g) and (h): calculated {(N2} — (Ni))/a 2 from upper panels. 



From equation l|17|) one can notice that if the second 
order N 2 in the exponential function can be neglected, 
ln(i?2i(^V)) is linear dependence on neutron number N 
for a fixed element with charge number Z, and a value 
can be expressed as (N 2 )/2a 2 2 — (N 1 )/2a 2 . Fig. |3]shows 
the centroid (N) and width a 2 of Gaussian isotopic distri- 
butions as a function of atomic number Z of the emitted 
products at different reaction time ( panel (a) - (f)) before 
and after secondary decay, panel (a) to (d), the difference 
of (N 2 ) in 48 Ca + 48 Ca system and (Aq) in 40 Ca + 40 Ca 
system increases when atomic number Z of the emitted 
products becomes large, in panel (e) and (f) of fig. |3 
the width a 2 of isotopic distributions does not same in 
two similar reactions too, it increases with the increas- 
ing of fragment atomic number Z ', reaches saturation at 
some point. In fig. E^g), ((N 2 ) - (N 1 ))/2a 2 is plotted as 
a function of Z, here a 2 = (a 2 + a 2 )/2 was adopted, it 
redisplays the a behavior in fig. its value raises from 
the light fragments to heavy fragment, especially in the 
heavy fragment part, it raises sharply. It decreases in 
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FIG. 4: The ratio of isoscaling parameters a between final 
and primary products. 



the region of light and intermediate mass fragments with 
increasing of the reaction time, but opposite trend in the 
region of heavy fragments. This reveals that neglecting 
the second order in equation (|17(l is reasonable in most 
cases. 

The secondary decay effect on isoscaling was investi- 
gated by the GEMINI code [H HU calculation. In this 
code, the fragment charge number Z, fragment mass A 
and the excitation energy of the fragment is the only in- 
put for the GEMINI code. In each event simulated by 
the IQMD, each fragment with charge number great than 
7 will consider sequential decay, the reason why the se- 
quential decay for fragment charge number great than 7 
will be discussed later. Since in the IQMD model, the 
angular momentum conservation was not taken into ac- 
count, the angular momentum of each fragment is not 
calculated, the angular momentum of the hot product 
was set to zero as input to GEMINI. Isoscaling behav- 
ior of final products after the secondary decay are also 
observed, a parameters were extracted with equation 
(|14f> as a function of the fragment atomic number Z and 
shown in fig. |2Ib), one can find that the dynamical ef- 
fect due to different reaction times is still clear for the 
light fragments, here the different reaction time refers to 
the time the fragments were collected in IQMD calcula- 
tion and the GEMINI calculation starts. The secondary 
decay effect, however, does not depend on the switch- 
ing time of dynamical IQMD to GEMINI code for the 
final products after secondary decay, a basically keeps 
not change for the light fragments even though the sec- 
ondary decay, but increases for the heavy products with 
the increasing of the products atomic number Z. In order 
to quantitatively compare the secondary decay effect on 
isoscaling parameters, the ratio between cxgemini and 
diQMD was plotted in fig. |3 one can see that the ratios 
fluctuate around 1.0 for the light fragment with atomic 
number Z < 8 and ~ 1.2 for the products charge number 
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FIG. 5: Comparison of charge distribution between the ex- 
perimental data (solid square) |3fj | and IQMD calculation 
(line) or IQMD + GEMINI (histogram) calculation at dif- 
ferent reaction time. 

Z > 8. 

The light fragments include both the primary prod- 
ucts which do not experience sequential decay and the 
secondary decayed products from the heavy primary 
products, however, this mixture of light fragments do 
not affect the isoscaling phenomenon. It reveals that 
the light primary fragments have same N/Z as the fi- 
nal products decayed from the heavy primary prod- 
ucts, they have achieved local equilibration or stable 
isotopic distribution on the isospin degree of freedom, 
so N/Z of light fragments is not affected by the sec- 
ondary decay process. It is a good observable to study 
the primary products isospin information, which have 
been confirmed b y many experiment and theoretical 
researches @ ©, GH 111 GS 13 111 El El El . 

In IQMD and GEMINI calculation, the fragment 
charge distributions of IQMD and GEMINI with the ex- 
perimental data |36| have been compared in fig. [5] The 
top panel is the comparison of fragment charge distribu- 
tion between experimental data [3(j and the IQMD cal- 
culation, it looks that heavier products have been overes- 



TABLE II: Average excitation energies per nucleon and the 
temperature calculated at different reaction times for two sim- 
ilar reaction systems 40 Ca + 40 Ca and 48 Ca + 48 Ca with in- 
cident energy Ei nc = 35 MeV/A and impact parameter 6 = 
lfm. 



t 


160/m/c 


200/m/c 


240/m/c 


320/m/c 


400/m/c 


(E*/A) a 


4.812 


4.524 


4.274 


3.871 


3.584 


rj-iO, 


6.937 


6.726 


6.537 


6.222 


5.987 


(E*/A) b 

rj-fb 


5.511 


5.214 


4.943 


4.496 


4.179 


7.423 


7.221 


7.030 


6.705 


6.464 


(TV 


7.180 


6.974 


6.784 


6.463 


6.226 



"average excitation energy per nucleon and temperature in colli- 
sion 40 Ca + 40 Ca 

''average excitation energy per nucleon and temperature in colli- 
sion 48 Ca + 48 Ca 

'average temperature (T) of 40 Ca + 40 Ca and 48 Ca + 48 Ca re- 
actions, (T) = (T a + T*)/2 

timated, but the intermediate mass fragments were un- 
derestimated in calculations. This problem was nicely 
solved after the secondary decay with the GEMINI code 
was taken into account (see bottom panel of fig. [5J, espe- 
cially when the switching time of IQMD fragments equals 
to 240 fm/c or 320 fm/c. The best description of the final 
charge distribution of products after GEMINI calculation 
lies on the bump of Z ss 6, this bump emerges only when 
secondary decay of the primary fragments with Z > 8 
are taken into account and while the secondary decay of 
primary fragments with Z < 8 are not considered. This 
may be the condition that the isoscaling parameters a is 
not affected by the secondary decay. 



IV. DENSITY DEPENDENCE OF THE 
SYMMETRY ENERGY COEFFICIENT 

The excitation energy per nucleon (E* /A) was calcu- 
lated from equation 1|13|) . E*/A distributions have the 
Gaussian shape, from which the average excitation en- 
ergy per nucleon ((E* /A)) can be obtained. In table UTI 
list the average excitation energy per nucleon at differ- 
ent reaction times for both similar reaction systems, the 
temperature is extracted from equation E*/A = a ■ T 2 , 
a = 1/10 was adopted. The excitation energy of the 
system has high values at early reaction time and it de- 
creases with the evolution of the system. At different 
reaction time the temperature extracted from two reac- 
tion systems 40 Ca + 40 Ca and 48 Ca + 48 Ca has approx- 
imately same values, so the average temperature in this 
two reactions are calculated simply to use as a parameter 
in following calculation, which are listed in table UTI too. 

Using equation (J2J, from the extracted values of a, av- 
erage temperature (T) of the two systems, and the {Z/A) 
of the fragments in two reaction systems, symmetry en- 
ergy coefficient C sym in the EOS can be derived. In 
present calculation, the (Z/A) was fitted from its Z/A 
distributions of the primary and final products respec- 



7 



28 




Atomic Number Z Atomic Number Z 



FIG. 6: Symmetry energy coefficients C sym as a function of 
the primary products atomic number Z from IQMD calcula- 
tion. 



FIG. 8: Average central density {pc/po) of the clusters in 
the rest frame of cluster as a function of the cluster proton 
number Z at different times for two reactions. 
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FIG. 7: Symmetry energy coefficients C sym as a function of 
the final products atomic number Z from IQMD+GEMINI 
calculation. 



tively, which are Gaussian distribution too and can get 
the average Z/A in two reactions {Z\/Ai) and (Z 2 /A 2 ). 
The derived symmetry energy coefficient C sym was plot- 
ted as a function of the atomic number Z of products for 
primary products in fig. Eland final products in fig. [7| 

One can find that the symmetry energy coefficient 
C S ym for different products with different atomic number 
Z is not a constant in the present calculation, it increases 
with the increasing of the products charge number Z, and 
the dynamics effect is almost vanished. To understand 
why the symmetry energy coefficient C sym was not a con- 
stant for different products, the densities of the products 
with different atomic number Z in IQMD simulation were 
extracted. In the IQMD frame, the fragments are treated 
as cluster by simple coalescence model, the central densi- 



ties of clusters were calculated by equation (@J with r = 
in the cluster coordinate phase space, and sum over all 
nucleons inside the same cluster. The central density dis- 
tribution pc/po of the clusters is generally Gaussian form 
and (pc/po) can be fitted from this Gaussian distribution, 
here po — 0.16/to~ 3 refers to the normal nuclear matter 
density. Fig. shows the fitted average central densities 
as a function of the cluster charge number Z at different 
reaction times in two reactions 40 Ca + 40 Ca and 48 Ca + 
48 Ca. The light clusters are always in the very low den- 
sity region, while the heavy clusters can approach to the 
saturation density. Remind that the nucleon can be con- 
sidered to be emitted if its density is lower than a certain 
low values, eg. po/10 in some transport calculations, the 
above curve of p(Z) support this argument. For reaction 
system 48 Ca + 48 Ca, a little higher average central den- 
sities (pc/po) of clusters than those in 40 Ca + 40 Ca are 
observed since the cluster with same proton numbers Z 
formed in reaction system 48 Ca + 48 Ca has more neu- 
trons than in reaction system 40 Ca + 40 Ca, which has 
been verified by the (N) distribution from fig. [21 The 
average central density of clusters in two reaction sys- 
tems have some divergence, in the later discussions, the 
average values of (p c /po)ave = ((p c /po)i + {Pc/pa)2)/^ of 
these two reactions is adopted. 

Since both symmetry energy coefficient C sym and the 
average cluster central density (p c / po) ave are correlated 
with the atomic number Z of products, the symmetry en- 
ergy coefficient C sym can relate with the average central 
density {p c / Po)ave by eliminating the products atomic 
number Z. The slight difference of different reaction 
time information was expressed by adding error bars to 
the derived C sym . Fig. ^3 demonstrates the relationship 
between symmetry energy coefficient C sym and average 
central density (p c /po)ave- The density dependence of 
the symmetry energy term in equation of state (EOS) 
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FIG. 9: The relationship between symmetry energy co- 
efficient C S ym and the average central density {p c / po)ave, 
solid square represents the IQMD calculations, the error bar 
represents the average over different reaction times. Panel 
(a): different lines represent the formula 11811 fitting with 
Co = 16.5 MeV and different parameter 7 which are shown 
in the insert; Panel (b): Combination of two set parameters 
C sym ({p c / P0 }) = 12 ■ (pc/po) ' 7 + 5 ■ (pc/po) 2 - 



FIG. 10: The relationship between symmetry energy coeffi- 
cient Csym and the average central density (p c /po)n«e, solid 
square represents the IQMD+GEMINI calculations and the 
error bar represents the average over different reaction times. 
Panel (a): different lines represent the formula l|18[l fitting 
with Co = 26MeV and different parameter 7 which are dis- 
played in the insert; Panel (b): Combination of two set pa- 
rameters C sym ({Pc/po)) = 12 ■ (Pc/po)"' 7 + 5 ■ (p c /po) 2 '° 



was generally expressed by following equation 

C sym (p) = C (— V (18) 

where 7 represents the stiffness of the symmetry energy. 
Some tests have been done for the determination of the 
parameter Co and 7 in fig. I1UI In panel (a) a constant Co 
= 16.5 MeV is used but some different values of 7, namely 
0.5, 0.7, 1.0, 1.3, 1.6 and 2.0 were taken, respectively, to 
fit the IQMD calculation results, one can find that 7 = 
1.0 seems to be the best fit to the IQMD calculation 
than other 7 values in the low (p c /po)ave < 1-0 region, 
but in (p c /po)ave > 1-0 region divergence arises between 
the IQMD calculation and the fitting, it seems that in 
high (p c /po)ave > 1-0 part, C sym raises sharply which 
corresponds to a 7 value great than 1.0. In panel (b) the 
line is two different parameters combination of equation 

(EH), C sym ({p c /p }ave) = 12 • (Pc/P0)al + 5 • (Pc/ Po)ave, 

in the low density region {p c /po)ave _ 1-0, has the 



similar shape with the above mentioned simple form 
C sym {(pc/po)ave) = 16.5 • (p c /po)L°e > but better behavior 
than the simple form with only 7 = 1.0. This demon- 
strates that in the whole density region, the symmetry 
energy coefficient C sym as a function of the nuclear den- 
sity can not be well described by only one simple equation 
(|18|l in the present model calculation, but expressed by 
the various contribution from two sets of EOS (|18|l . In 
present calculation of IQMD, the simple equation such 

as C sym ((p c /p ) ave ) = 12 • (pc/Po)°ave + 5 • (Pc/Po) 2 ave & a 

good approximation for the C sym and (p c / Po}ave- 

Fig. plots the relationship between C sym and 
(pc/po)ave after considering the secondary decay. Panel 
(a) shows the IQMD + GEMINI calculation C sym as 
a function of average central density {p c J Po)ave, several 
fits with a constant Co = 26MeV and different 7 = 0.5, 
0.7, 1.0, 1.3, 1.6, 2.0 of equation (|18|l have been made, 
in which 7 = 1.3 looks to fit the calculation best in 
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the density region (p c /po)ave < 1.0, equation (|18(l with 
7 >1.3 works better in the higher density part. Equation 

C sym {(pc/Po)ave) = 16 • (Pc/Po)ave + 13 " (Pc/Po)lve wa S 

also adopted in panel (b) to fit the calculation in fig. 1101 
it can give a good coincidence between the calculation of 
IQMD+GEMINI and the equation (HHJl fitting. 

No matter a simple form of equation i|18fl was used, 
or combination of different set parameters of equation 
(|18fl was used, different parameter Co was got for the 
primary products in IQMD calculation and final products 
in IQMD+GEMINI calculation. In the simple form of 
equation l|18|) for fitting IQMD calculation, Co = 16.5 
MeV was got, while fitting IQMD + GEMINI calculation 
Co = 26 MeV was got. The different symmetry energy 
coefficients Co extracted from primary and final products 
indicate that the secondary decay affects the extraction 
of symmetry energy coefficient Co- 

From the fits to the IQMD and IQMD+GEMINI cal- 
culations, the subnormal density region can be described 
better by a soft symmetry energy with 7 < 1.0, the su- 
pernormal density region can be described better by a 
stiff symmetry energy with 7 > 1.0. This demonstrates 
that in the subnormal and supernormal density regions, 
the symmetry energy shows different behavior, the best 
description of the relationship between symmetry energy 
coefficient C sym and the nuclear density (p/po) in full 
density region should be combination of equation 118|) 
with at least one "soft" symmetry term with 7 < 1.0 and 
one "stiff" symmetry term with 7 > 1.0 28]. 

V. CONCLUSIONS 

In summary, we applied the Isospin dependent Quan- 
tum Molecular Dynamics (IQMD) model to investigate 
the isoscaling behavior in dynamical and statistical se- 
quential decay processes of both reactions 40 Ca + 40 Ca 
and 48 Ca + 48 Ca at b = 1 fm. The calculation illus- 



trates that both primary and final isotopic yields show 
the isoscaling behavior. 

The isoscaling parameter a have similar values for the 
light and intermediate mass fragments, but it increases 
much with the fragments become heavy, a shows some 
dynamical effect, i.e. it drops slowly with the increasing 
of reaction times. The light final products which do not 
experience secondary decay itself are not affected by the 
secondary decay, but the isoscaling parameter a of the 
heavy products which experience secondary decay pro- 
cess increase in comparison with the primary products. 

Average central density of the clusters is calculated 
in the IQMD calculation, it raises with the increasing 
of charge number of clusters until the saturation den- 
sity close to the normal nuclear matter density po = 0.16 
fm -3 . From the equation (|18|) . the density dependence of 
the symmetry energy coefficient C sym can be extracted 
from 'the IQMD and IQMD+GEMINI calculation. The 
results show that even though the density dependence 
of the symmetry energy coefficient C sym can be approx- 
imately described by a 7 ~ 1.0, it is more reasonable 
to express this dependence by the combination of " soft" 
and "stiff" symmetry energy term. From equation i|18fl . 
the extracted symmetry energy coefficient C syrn is also 
affected by the secondary decay process, i.e. it leads to 
different values of Co and 7 from the initial values of 
IQMD. 
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